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MODELING  OF  III-V  FIELD  EFFECT  TRANSISTORS 


Abst  ract 


A two-dimensional  time  domain  computer  program  has  been 

formulated  for  investigating  device  physics  of  III-V  field 

effect  transistors.  The  program,  with  our  modified  version 

( 1 ) 

of  Hockney's  algorithm  for  solving  the  two-dimensional 
finite  difference  Poisson's  equation,  is  about  two  orders  of 
magnitude  faster  than  comparable  programs  using  over-relaxation 
techniques^  ^ . The  increase  in  computational  speed  makes  it 
practical  to  simulate  time  domain  behavior  of  the  device  in 
details.  Internal  mobile  carriers  and  electrostatic  potential 
distributions  and  the  drain  current  can  be  examined  as  functions 
of  time  when  the  device  is  subjected  to  a time-step  input  signal. 
Relevant  information  is  then  derived  from  the  large  signal  impulse 
response . 


Using  this  technique  we  investigated  the  effect  of  gate 
length  to  channel  width  ratio,  high  field  domain  formation, 
tailoring  of  doping  profile,  vertical  vs  planar  FETs  and 
comparison  of  the  GaAs  and  InP  FETs.  Additionally,  the  device 
cut-off  frequency  and  intrinsic  delay  time  as  functions  of  gate 
length  were  determined.  Significant  findings  are:  1)  The 
dynamic  transfer  characteristic  (drain  current  vs.  time-step 
gate  voltage)  of  the  FET  is  like  that  of  a series  R-C  circuit. 
2)  The  device  behavior  depends  strongly  on  the  gate  length  to 
channel  depth  or  L^/d  ratio.  3)  For  planar  devices  with  a 
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Lq/cI  < 4 the  static  ID-V[)  characteristics  have  a pronounced  slope 
beyond  current  saturation  due  to  carrier  injection  into  the  semi- 
insulating  substrate.  4)  The  fj  and  the  device  intrinsic  delay 
time  are  not  determined  by  the  metallurgical  gate  length  alone 
due  to  the  two-dimensional  effect  creating  depletion  regions 
extending  beyond  the  gate.  The  simulation  of  InP  FET  was  carried 
out  over  a limited  range  of  operation. 


Objectives  and  Accomplishments 


The  objectives  of  the  past  year's  efforts  are  the  following: 
l)  Develop  the  model  and  formulate  the  relevant  mathematical 
equations  needed  to  obtain  a quantitative  description  of  field- 
effect  transistor  physics,  2)  Modify  the  existing  subroutines 
developed  under  contract  F44620-76-C-0043  in  order  to  include 
time-dependent  carrier  response  to  electric  field,  carrier  trapping 
and  more  flexible  boundary  conditions  for  the  simulation  of 
different  device  configurations,  3)  Simulate  and  study  device 
physics  of  planar  GaAs  and  InP  FETs  and  4)  Simulate  and  study 
device  physics  of  vertical  GaAs  and  InP  FETs. 


All  the  objectives  with  the  exception  of  the  study  of 
effects  of  carrier  trapping  have  been  fully  achieved.  In  the 
sections  that  follow  we  describe  (a)  the  model,  the  physical 
mechanisms  and  the  mathematical  equations  used  in  the  computer 
program,  (b)  modification  of  Hockney's  algorithm  and  the  sub- 
routine developed  under  F44620-76-C-0043  for  simulations  of 
planar  and  vertical  FETs  (c)  major  results  of  our  investigations. 
A brief  summary  is  given  in  section  (d) 
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( a ) Ihe.Comgut er_P rogram 

The  physical  models  for  the  planar  (one-sided)  and  vertical*  * 

(two-sided)  are  shown  in  Figure  1.  The  device  internal  physical 

I, 

effects  treated  in  our  program  are  the  following:  two-dimensional 
carrier  transport  from  electric  field  and  diffusion;  carrier  and 
charge  conservation  during  transport  and  avalanche  process;  space 
charge  effects  resulting  from  mobile  and  fixed  charges;  field  and 
time  dependent  avalanche  process;  nonlinear  field-dependent 
carrier  velocities  and  velocity  over-shoot  when  applicable.  The 
external  boundaries  are:  ohmic  contacts  on  the  source  and  drain, 

Schottky  barrier  on  the  gate  electrode,  electrode  voltages,  and 
conservation  of  terminal  currents. 


The  equations  governing  the  two  dimensional  device  are: 
Poisson's  equation 

V24>  2 - - ( Nd  + p - n ) (1) 


being  simulated.  However,  instead  of  the  electric  field  alone, 

*We  use  the  term  'vertical'  to  suggest  the  actual  construction 
of  such  dual  gate  F.E.T.s  on  a substrate  which  would  serve  as 
the  source. 


ki 


-5- 


the  force  field  is  defined  to  include  diffusion  so  that 


for  electrons 


and 

T 

f = - V 4*  ~ — VD  for  holes. 

P P 

These  forces  as  defined  are  consistent  with  the  concept  of  the 
quasi-fermi  potential  and  are  used  in  place  of  the  electric 
field  in  the  v-E  characteristics.  In  this  way  the  carrier 
saturation  velocity  cannot  be  exceeded,  which  would  be  the  case  if 
the  electric  and  diffusion  forces  are  treated  separately. 

To  obtain  good  resolutions  in  the  simulation,  the  two- 
dimensional  device  is  devided  into  96x48  cells  as  shown  in 
Figure  2.  Equations  (1)  and  (2),  in  finite  difference  form,  are 
computed  using  the  explicit  technique.  Stable  solutions  of 
the  transport  equation  are  ensured  by  choosing  the  time  increment 
(At)t  (subscript  t for  "transport")  such  that 


( At ) t 


v < 
max  - 


1 


(3) 


Ax  T 

where  v is  the  maximum  velocity  a charge  carrier  can  achieve. 

ID  3 X 

For  computer  economy,  a larger  time  increment  (At)p  =m  (At>t  is 
used  for  updating  the  Poisson's  equation.  The  value  of  the 
integer,  m , is  determined  by  the  size  of  C A t ) t and  the  device 
time  constant.  Typically,  m ranges  from  2 to  8. 
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(b)  Modi^f2cat2on_of_Prey20us_Sub  — routines 

Sub-routines  developed  under  AFOSR  contract  F44620-76-C-0043 
for  bipolar  device  simulation  were  modified  for  FET  applications. 
However,  one  of  the  sub-routines,  the  Hockney  fast  algorithm  for 
solving  Poisson's  equation,  was  found  not  to  be  compatible  with 
the  boundary  conditions  required  to  describe  the  FET.  This  dif- 
ficulty was  resolved  by  incorporating  an  innovation  that  is  dis- 
cussed here. 

The  FET,  as  shown  in  Figure  2,  contains  both  Dirichlet  and 

Neumann  boundaries,  while  the  Hockney  technique  is  applicable 

to  either,  but  not  mixed,  boundaries.  In  addition,  the  extent 

of  the  boundary  is  time  dependent.  To  alleviate  the  difficulties, 

the  Newmann  boundary,  which  arises  from  the  gate  depletion  region, 

is  converted  into  Dirichlet  in  each  time  increment  (At)p.  This 

could  be  done  because  no  current  can  flow  across  the  exterior 

boundary  (x=a)  of  the  depletion  region;  thus  the  condition  4^1  = 0 

®x  * x = a 

prevails.  This  condition,  together  with  the  known  doping  density 
and  potentials  at  all  electrodes  permits  the  evaluation  of 
4>l(a*y*t)  and  <t>2^a,y,t>  in  Figure  2 via  the  one-dimensional  Poisson 
equation  in  y.  Thus,  the  potential  $(a,y,t)  for  a <_  y < b is 
completely  specified.  The  function  4>(a,y,t)  is  Fourier-analyzed 
for  each  (At)p.  The  resultant  Fourier  components  are  then  applied 
to  appropriate  locations  in  Hockney's  algorithm.  This  procedure 
is  repeated  for  the  boundary  at  x=0  for  the  vertical  FET.  The 
boundaries  at  y=0  and  y=b  are  considered  as  planes  of  even 
symmetry  appropriate  to  the  cosine  expansion  used  in  the  Fourier 
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analysis.  Such  symmetry  property  implies  that  the  Ey  is  zero 
at  these  edges  and  does  not  produce  error  as  long  as  the  source 
and  drain  electrodes  extend  sufficiently  inward. 

(c)  Ma  j_o  r_ R e s u 1 1 s 

At  the  inception  of  present  work,  it  was  decided  that 
instead  of  a sinusoidal  signal,  a time  step  input  would  be  used 
in  probing  the  device.  As  will  be  seen  later,  the  large  signal 
impulse  response  produced  more  information  on  the  device  physics. 
Furthermore,  a time  step  signal  is  consistent  with  logic  device 
operation  of  the  FET,  which  is  of  major  interest  at  present. 

In  the  simulation,  the  device  was  connected  in  the  common 
source  configuration  with  a fixed  drain  bias.  The  gate  voltage 
VgS  was  positive  and  equal  to  the  Schottky  barrier  potential  for 
t<  0.  At  t=0,  the  gate  voltage  was  switched  instantaneously 
to  a negative  value.  The  mobile  carrier  and  electrostatic 
potential  distributions  and  the  drain  particle  current  were 
then  recorded  as  functions  of  time  starting  at  t=0.  For  computer 
economy,  a doping  density  of  5 x 10^  cm“3  in  the  active  region 
and  a channel  depth  of  1yM  were  used.  For  our  purposes,  these 
parameters  do  not  restrict  the  validity  of  the  computational 
results;  the  result  for  higher  doping  level  can  be  scaled  from 
that  for  the  5 x lO1^  doping.  The  main  physical  variable  in 
the  investigation  was  the  gate  length.  The  major  results  are 
presented  in  the  following  sub-sections. 
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(1)  Imgulse_Resgonse 

The  drain  currents  under  impulse  conditions,  for 
three  different  gate  lengths,  are  shown  in  Figure  3.  The 
devices  simulated  had  a pinch-off  voltage,  Vp  of  3.77V, 
according  to  the  gradual  channel  approximation,  while  the 
gate  input  was  -4 V.  In  the  long  gate  case  (Lg/d=8),  ex- 
cept for  the  slight  Gunn  oscillation,  the  semi-log  plot 
of  Io  versus  t is  indeed  a straight  line.  Thus,  the  dy- 
namic transfer  characteristic  of  the  FET  can  be  repre- 
sented by  a series  R-C  circuit.  The  device  time  constant 
t determined  from  tf>e  slope  is  related  to  the  transcon- 
ductance cut-off  frequency  by  the  relation^^ 


Effects  of  Lg/d  ratio  are  apparent  when  the  drain 
current  waveform  is  inspected.  We  shall  discuss  the  gate 
length  mostly  in  terms  of  this  ratio  rather  than  the  ab- 
solute value  of  Lg.  The  dividing  line  between  a long 
gate  and  a short  gate  device  appears  to  be  Lg/d=4.  The 
short  gate  devices  have  two  distinct  features:  softer 
pinch  off  at  Vg^Vp  and  larger  degree  of  current  oscilla- 
tion. The  time  constant  of  short  gate  devices  was  de- 
termined from  the  initial  slope  of  the  drain  current  as 
indicated  in  Figure  3.  For  large  values  of  t,  the  drain 
current  becomes  saturated  permitting  the  study  of  steady 


state  behavior 
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A critical  observer  may  note  that  the  magnitude  of 
the  initial  current  in  Figure  3 does  not  scale  inversely 
with  the  length.  This  could  be  accounted  for  by  the 
geometrical  effect  on  the  electric  field  and  the  peculiar 
velocity-field  characteristic  of  GaAs.  In  the  following 
sub-sections  the  dynamic  behaviors,  drived  from  impulse 
response,  are  described. 

(2)  Gunn_Os c i_  1 1 a t i on 

Formation  of  a high  field  domain  in  the  gate-drain 
region  has  been  observed  in  two-dimensional  computations 
C5-9D.  Grubin^)  pointed  out  that  the  high  field  domain 
is  not  necessarily  of  negative  differential  mobility  in 
origin,  rather  it  could  also  arise  from  the  geometrical 
effect  within  the  FET.  Additionally,  Gunn  oscillation 
has  also  been  reported.  Grubin^^  found  that  the  oscilla- 
tion can  be  suppressed  by  increasing  the  drain  bias  voltage 
while  Yamaguchi  et  al.^9)  defined  a stable  region  where  the 
channel  depth  is  reduced  to  certain  values.  Our  computer 
simulation  not  only  supports  these  findings  but  also 
illustrates  how  and  where  the  Gunn  region  exists. 

Since  the  high  field  domain  is  of  interest  in  under- 
standing the  GaAs  FET  device  physi  c s ^ ^,1  1)  , we  show  our 


computational  results  in  some  details  here.  The  time 
evolution  of  the  mobile  charge  distributions  in  a short 
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gate  (Lg/d=2)  device  under  impulse  conditions  are  shown 
in  Figure  4.  Although  hundreds  of  carrier  distribution 
frames  are  computed  for  one  computer  run,  only  pertinent 
snap  shots  are  displayed  to  illustrate  the  device  behavior. 
The  undulation  of  the  conducting  channel,  as  shown  in  frames 
from  t=40  picoseconds  to  t=68  picoseconds,  is  two-dimensional 
and  Gunn  effect  in  origin.  These  frames  cover  a time 
duration  in  the  current  saturation  region  (See  Figure  3) 
and  thus  the  oscillation  is  a steady  state  as  well  as  tran- 
sient behavior.  At  t = 4 0 psec.  the  channel  is  widest  near 
the  drain.  As  the  bulge  disappears  into  the  drain,  a new 
one  is  initiated  near  the  source.  The  cycle  is  complete 
at  t=68  psec.  The  Gunn  oscillation  is  further  demonstrated 
by  the  potential  distribution  in  Figure  5.  By  examining 
the  time  dependence  of  the  equi pot ent i a l lines,  it  is 
apparent  that  the  most  active  region  is  in  the  conduction 
channel.  This  accounts  for  the  strong  Gunn  oscillation  in 
short-gate  device  since  the  drain  voltage  penetrates  deeply 
into  the  channel.  The  steady  state  potential  distribution 
in  a long-gate  device  is  also  given  in  Figure  5 for  compari- 
son. The  main  body  of  the  high  field  region  appears  to  be 
"stationary"  except  for  the  slight  shifting  of  the  equipo- 
tentials. 


Increasing  the  drain  voltage  diminishes  but  does  not 
completely  suppress  the  Gunn  oscillation,  as  demonstrated 
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by  the  current  waveform  in  Figure  6.  Such  undulating 
field  region  is  not  determined  by  the  absolute  value  of 
the  channel  depth  alone,  rather  it  is  a function  of  Lg/d 
ratio  as  illustrated  in  Figure  3.  This  is  consistent 
with  Yamaguchi's  results  since  his  stable  region  is 
found  at  Lg/d>4.  Although  Gunn  oscillation  can  be  mini- 
mized (5,9),in  all  our  calculations  it  never  completely 
disappears  even  after  a very  long  time  (t>10x). 

The  fact  that  the  drain  current  may  appear  to  be  relatively 
quiescent  while  the  internal  carriers  go  through  two-dimensional 
undulation  is  because  the  drain  current  is  an  integral  of 
the  carrier  motions.  We  suggest  that  the  so-called 
"stationary"  domain,  based  on  external  current  observation, 
may  only  be  an  approximation  to  what  actually  transpires 
internally  in  the  GaAs  FET.  It  is  noted,  however,  that 
our  computations  are  far  from  exhaustive. 

(3)  Cut-0ff_Freguency_and_Sw2tch2ng_Dela£ 

Numerical  result  of  the  transconductance  cut-off 
frequency  as  a function  of  gate  length  for  the  GaAs  FET 
is  shown  in  Figure  7.  The  values  were  determined  from 
the  impulse  response  drain  current  and  equation  (4). 

Carrier  velociy  over-shoot  was  neglected  in  these  calcu- 
lations because  it  was  found^0,11)  that  this  effect  is 


I l 


insignificant  for  gate  length  of  1.0  ym  and  longer. 
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In  an  actual  device,  carriers  are  injected  into  the  source 
region  through  an  ohmic  contact  and  traverse  a finite  source 
region  before  reaching  the  gate.  Thus,  the  overshoot  effect 
is  not  as  significant  as  calculated  in  the  case  of  one- 
dimensional approximation  and  using  a more  abrupt  electric 
field  slope. 

For  short-gate  cases,  velocity  saturation  takes  place 
in  the  conducting  channel  and  f0~fy,  the  current  gain  cut-off 
frequency.  Although  for  a given  device,  fj  depends  on 
operating  pa ramete rs ^ 2) , our  results  compare  favorably 
with  published  experimental  data^^-15).  Thus,  fy  as 
determined  from  the  impulse  response  corresponds  to  that 
deduced  from  measurements  at  V g 5 = 0 . It  should  also  be 
noted  that  the  two-dimensional  simulator  accounts  for 
depletion  regions  at  both  ends  of  the  gate.  Consequently, 
the  effective  gate  length  is  appreciably  longer  than  the 
metallurgical  value  as  suggested  by  Engelmann  and  Liechti^^). 
The  effective  gate  length  probably  accounts  for  the  differ- 
ence between  our  fy  and  the  one-dimensional  results  of 
Maloney  and  Frey^^). 

The  intrinsic  device  switching  delay  was  also  deter- 
mined from  the  i p-t  characteristics.  The  delay  time  is 
defined  as  the  duration  in  which  the  drain  current  is 
reduced  to  10%  of  its  intitial  value.  For  nonpinched  off 
cases  such  as  in  Figure  3 the  initial  slope  was  extrapolated 
to  the  10%  level.  The  delay  time  vs.  gate  length  is 
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* shown  in  Figure  8.  Unfortunately,  there  is  only  one 
experimental  data  point^-^  to  compare  with  our  results. 
According  to  Figure  8,  the  GaAs  FET  appears  to  be 
intrinsically  a very  fast  logic  device. 

(4)  §ubst rat e_Cur rent 

As  demonstrated  earlier,  short  gate  (Lg/d<4)  GaAs 
FETs  are  prone  to  Gunn  oscillation  and  require  a larger 
gate  voltage  to  effect  pinch-off.  The  mechanism  for 
these  peculiarities  is  the  deep  penetration  of  drain 
voltage  into  the  conducting  channel  as  shown  in  Figure  5. 

• In  planar  short  gate  FETs  the  deep  penetration  also  causes 
carriers  in  the  conducting  channel  to  be  injected  into 

the  sem i - i n su l a t i ng  substrate.  The  effect  is  demonstrated 

in  Figure  9 where  carrier  densities  in  the  substrate  of  a 

short-  and  a long-gate  devices  are  compared.  It  is  apparent 

that  the  substrate  current  is  significant  in  the  short-gate 

but  not  so  in  the  long-gate  device.  Since  the  origin  of 

carrier  injection  is  in  the  drain  voltage,  we  suggest  that 

this  is  the  mechanism  for  the  pronounced  rising  slope  of 

the  drain  current  beyond  the  saturation  point  in  sub-micron 

gate  (L  <1/4  ym)  devices^^.  As  the  gate  length  is  made 

n 

smaller,  it  becomes  increasingly  difficult  to  maintain 
uniform  Lq  and/or  d in  practice  and  hence  the  soft-pinch 
off  behavior  is  determined  by  a region  with  the  lowest  Lg/d 
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ratio.  Thus,  our  result  implies  that  there  might  be 
a practical  lower  limit  on  the  gate  length  base  on  device 
performance  with  current  technology.  Similar  computation 
results  on  substrate  current  have  been  reported  by  Reiser^®). 
However,  he  did  not  consider  it  as  a function  of  gate 
length. 

(5)  lit 

Tailoring  of  doping  profile  along  the  gate  width  has 
been  reported  by  Siliger  and  Ward^9).  With  a certain 
variation  of  N&  as  a function  of  z (direction  along  the 
gate  width)  the  linear  dynamic  range  of  the  device  was 
increased  by  several  orders  of  magnitude. 

We  found  that  tailoring  of  doping  profile  along  the 
gate  length  (y  direction  in  Figure  2)  can  significantly 
increase  the  cut-off  frequency  fQ.  A 2 pM-gate  device 
with  an  exponentially  tapered  profile  has  an  f 0 = 1 3 . 5 GHz 
as  compared  to  9 GHz  for  the  uniform  profile.  In  this 
device  the  doping  density  varied  from  10^  cm~3  a*  the 
source  to  10^  Cm_^  at  the  drain.  The  increase  in  f0  is 
due  not  only  to  the  higher  conductivity  in  the  source 
region,  but  also  to  the  net  lowering  of  the  o ite-to-source 
distributive  capacitance. 

A one-step  doping  variation  along  the  gate  produced 
an  even  more  dramatic  result.  The  doping  level  is 
1 x 1D^cm-3  from  source  to  the  middle  of  the  gate  and 
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1 x 10^cm“3  from  there  to  the  drain.  This  doping  scheme 
effectively  moves  the  source  to  the  middle  of  the  gate, 
thus  increasing  the  f0  to  17  GHz.  It  is  not  unreasonable 
that  the  f is  higher  than  the  1 ym-gate  device  (see 
Figure  7)  since  the  effective  channel  length  is  actually 
shorter.  A comparison  of  impulse  response  drain  current 
for  the  three  types  of  doping  profiles  is  given  in 
Figure  10. 

(6)  Ga As_vsi_InP_FETs 

Barrera  and  Archer^O)  found  that  for  a one-yM  gate 
device  the  InP  FET  has  an  fj  of  20  GHz  as  compared  to 
13  GHz  for  GaAs.  They  attribute  the  difference  to  the 
higher  peak  carrier  velocity  in  InP.  Turner  and  Wilson^1^ 
also  predict  that  cut-off  frequency  ratio  of  the  two 
materials  is  1.3  in  favor  of  InP. 

Our  computer  simulation  results  so  far  have  not  been 
adequate  to  confirm  the  above  findings.  The  impulse 
responses  of  GaAs  and  InP  FET  of  identical  construction 
(2  yM-gate)  and  operated  under  identical  conditions  are 
shown  in  Figure  11.  Similar  results  for  the  4 yM  gate, 
not  shown,  are  also  obtained.  The  notable  feature  is 
the  more  pronounced  oscillation  in  the  InP  drain  current; 
the  macroscopic  feature  of  current  decay  is  about  the  same. 
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In  actual  experiments  such  as  in  Ref.  (20),  operating 
conditions  can  readily  be  varied  to  achieve  best  device  per- 
formance. The  electric  field  at  which  the  velocity  peaks  is 
*10  kV/m  for  InP,  and  the  peak  velocity  is  larger  as  well  as 
broader  than  that  in  GaAs.  It  is,  therefore,  feasible  that 
the  carrier  traverse  a large  fraction  of  the  conduction 
channel  at  near  the  peak  velocity  when  the  device  parameters 
and  operating  conditions  are  appropriate.  Our  computer 
simulations  have  not  covered  sufficient  range  of  conditions 
to  permit  any  conclusion  on  the  fj.  Rather,  they  only 
indicate  that  the  InP  FET  is  more  prone  to  current  oscilla- 
tion. 

(d)  Summary 

We  have  described  the  formulation  of  a time-domain 
two-dimensional  computer  program  which  is  very  fast  in 
the  execution.  The  program  resolves  the  FET  into  96x48  cells 
and  simulates  its  dynamical  behavior  from  transient  to  steady 
state.  The  time  increment  used  for  up-dating  the  internal 
dynamical  quantities  is  small  compared  to  the  device  time 
constant  C(At)t  * . The  fine  spatial  and  temporal  resolu- 

tions permit  examination  of  the  FET  behavior  in  minute  details. 

Using  this  program  the  device  physics  was  studied  by 
observing  its  large  signal  impulse  response.  The  resultant 
internal  carrier  density  and  potential  distributions  indicate 
persistent  Gunn  oscillation  and  significant  substrate  current 
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for  short-gate  devices.  The  latter  explains  the  soft  pinch-off 
characteristics  of  sub-micron  gate  devices.  The  computed  fj 
and  intrinsic  switch  delay  time  agree  well  with  available 
experimental  data.  The  two-dimensional  effects  on  the  fT  was 
found  to  be  significant.  Finally,  tailoring  of  the  doping 
profile  can  increase  device  gain-bandwidth  product  at  the  expense 
of  its  power  capability. 


-18- 


Li_st_of_References 

1)  R.  W.  Hockney,  "A  Fast  Direct  Solution  of  Poisson's 
Equation  Using  Fourier  Analysis",  J.  Assoc,  Comput  . 
Mach.  12  PP.  95-113  (1965). 

2)  C.R.  Jershope,  "Numerical  Solution  to  2-dimensional 
Time-Dependent  Semiconductor  Equations",  Elect.  Lett 
20,  PP  431-433  (1975). 

3)  M.  Adler,  Private  Communication 

4)  J.R.  Hauser,  "Fundamentals  of  Silicon  Integrated 
Device  Technology",  (Eds.  R.M.  Burger  and 

R.P.  Donovan),  Volume  2,  Prentice  Hall,  1968. 

5)  H.L.  Grubin,  "Large  Signal  Numerical  Simulation  of 
Field  Effect  Transistor"  6th  Biennial  Conference 
on  Active  Microwave  Semiconductor  Devices  and 
Circuits,  Cornell  University,  1977. 

6)  D.P.  Kenedy  and  R.R.  O'Brien,  "Computer  Aided  Two- 
Dimensional  Analysis  of  the  Junction  Field  Effect 
Transistor",  IBM  J.  Res.  Develop.  14,  PP  95-116,  1970. 

7)  B.  Hinsworth,  "A  Two-Dimensional  Analysis  of  Gallium 
Arsenide  Junction  Field  Effect  T;ansistor  with  Long 
and  Short  Channels",  Solid  State  Electron.  15 

PP  1353-1361,  1972. 


8)  B.  Hinsworth,  "A  Two-Dimensional  Analysis  of 
Indium  Phosphide  Junction  Field  Effect  Transistor 
with  Long  and  Short  Channels",  Solid  State  Electron, 

16  PP  931  - 939,  1973. 

9)  K.  Yamaguchi,  S.  Asai,  and  H.  Kodera,  "Two-dimensional 
Numerical  Analysis  of  Stability  Criteria  of  GaAs 
FETs",  JEEE  Trans.  ED-23,  PP  1283-1290,  1976. 

10)  S.  P.  Yu  and  W.  Tantraporn,  "Time-Domain  Two- 
Dimensional  Computer  Simulation  of  T h re e -T e r m i na  l 
Semiconductor  Devices"  6th  Biennial  Conf.  on  Active 
Microwave  Semiconductor  Devices  and  Circuits,  Cornell 
University  1 977 . 

11)  M.  Shur,  "Influence  of  Nonuniform  Field  Distribution 
on  Frequency  Limits  of  GaAs  Field  Effect  Transistors" 
Elect.  Lett  12,  615  (1976). 

12)  R.W.H.  Engelmann  and  C.  A.  Liechti,  "Bias  Dependence 
of  GaAs  and  InP  MESFET  Parameters"  IEEE  Trans  ED-24, 

PP  1288-1296  (1977). 

13)  Van  Tuyl  and  Liechti,  "Gallium  Arsenide  Digital 
Integrated  Circuits",  Tech.  Report  AFAL-TR-70-40 
(1974). 

14)  H.  Q.  Tsung,  V.  Sokolov,  H.M.  Macksey  and  W.R.  Weisseman, 
"Microwave  Power  GaAs  FET  Amplifier"  IEEE  Trans. 

MTT-24 , PP936-943  (1  976). 

15)  L.  Napoli  and  R.  DeBrecht,  "Performance  and  Limitation 
of  FETs  as  Microwave  Power  Amplifier"  IEEE  Microwave 
Symp.,  Univ.  of  Colorado  (1973). 


"20* 


16)  T.  J.  Maloney  and  J.  Frey,  "Frequency  Limits  of 

GaAs  and  InP  Field  Effect  Transistors  at  300°K 
and  77°K  with  Typical  Active  Layer  Doping"  IEEE 
Trans.  ED-23,  PP  519  - (1976). 

17)  R.S.  Bullin,  et  al.,  "200  Nanometre  Gate  Length 
FETs"  6th  Biennial  Cornell  Active  Microwave  Semi- 
conductor Devices  and  Circuit  Conference,  Cornell 
Univ.  (1977). 

18)  M.  Reiser,  "Two-Dimensional  Analysis  of  Substrate 
Effects  in  Junction  FETs"  Elec.  Lett.  6,  PP  493-494 
(1970)  . 

19)  R.  L.  Seliger  and  J.  W.  Ward,  "Fabrication  of 
Lateral  Doping  ProfiLes  by  a Computer  Controlled 
Focused  Ion  Beam",  J.  Vac.  Sci.  Technol.  \2  PP.  1378  - 
1381  (1975). 

20)  J.  S.  Barrera  and  R.  J.  Archer,  "InP  Schottky  Gate 
Field  Effect  Transistor"  IEEE  Trans  ED-22  PP  1023- 
1030  (1975). 

21)  J.  A.  Turner  and  B.  L.  H.  Wilson,  "Implication  of 
Carrier  Velocity  Saturation  in  a Gallium  Arsenide 
Field  Effect  Transistor"  1968  Symp.  on  GaAs  Proceedings 
PP  195-204,  (Inst,  of  Phys.,  Phys.  Soc.  London,  1969). 


NORMALZIED  DRAIN  CURRENT 


Nd=5xI015  cm'3 
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Figure  3.  Semi-Log  Plot  of  Drain  Currents  Under  Large  Signal  Im 
pulse  Condition.  Note  that  the  short-gate  Lg/d  = 2 
device  is  not  pinch-off  and  has  a stronger  Gunn  oscil 
lation.  The  initial  slope  of  the  macroscopic  time  de- 
pendence (represented  by  dashed  lines)  is  used  in  the 
calculation  of  x . 
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SOURCE  GATE  DRAIN 


Figure  5.  a,  b,  and  c:  Electrostatic  Potentia.  Distributions  in 
the  LG/d  = 2 Device.  The  degree  of  penetration  of  the 
drain  voltage  into  the  conducting  channel  is  correlated 
to  magnitude  of  Gunn  oscillation  in  the  drain  current, 
as  well  as  the  time  variation  (indicated  in  ps)  of  the 
equipotential  lines.  Also  shown  (d)  is  the  potential 
distribution  of  a long-gate  device.  Note  the  scale 
distortion.  The  full  vertical  side  is  2pm  for  all 
frames  while  in  the  horizontal  direction  the  gate 
length  is  2pm  in  a,  b,  & c and  is  8pm  in  d. 
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Figure  6.  Impulse  Response  of  the  Drain 

Current  Showing  Effects  of  Drain 
Voltage  on  f0  and  Gunn  Oscillation 
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Figure  7 . Computed  Cut-off  Fre- 
quency vs.  Gate  Length: 
Present  Work  and  the 
Result  of  a One- 
dimensional Calculation 
of  Maloney  and  Frey. 
Experimental  data  points 
are  taken  from  the  ref- 
erences cited. 
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Figure  8.  Device  Time  Constant  T and  Intrinsic  Switching  Delay  Time  as  Function 
of  Gate  Lengtn.  One  experimental  point  from  van  Tuyl  and  Liechti  (12) 
is  shown. 
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(a)  (b) 

Figure  9.  Steady  State  Distribution  of  Mobile  Carriers  in  Short- 

and  Long-Gate  Devices.  (a)  Vertical  frame  = 5um,  active 
layer  thickness  d = 1pm  with  doping  of  5 x lO^cm""3  on 
4pm  insulating  layer.  Horizontal  frame  = 8pm  with  gate- 
length  Lg  = 2pm.  (Lg/d  = 2)  S,  G and  D denote  source,  gate 
and  drain  electrodes.  (b)  Same  as  (a)  except  horizontal 
frame  = 32pm  and  Lg  = 8pm.  Note  the  larger  degree  of 
carrier  injection  into  the  substrate  in  (a). 
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Figure  10.  Impulse  Response  of  the  Drain  Current  for  Three 
Identical  Devices  Except  the  Doping  Profile 
ND (y)  . (a)  Constant  5 x 10l5cm-3.  (b)  Exponen- 

tially decreasing  from  1 x 10^6  at  the  beginning 
to  1 x 1015  at  the  end  of  the  gate,  (c)  Abrupt 
change  from  1 x 10^  to  1 x 10*^  at  the  middle  of 
the  gate.  Other  device  parameters  are:  d = 1pm, 
Lg  = 2pm,  Vq5  = -4  volts,  V&S  = 10  volts. 
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Figure  11.  Impulse  Response  of  the  Drain  Current  for  Identical 

Devices  Except  the  Material  Velocity-field  Character- 
istics. Device  parameters  and  the  operation  conditions 
are  the  same  as  in  Figure  10a. 
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